High-Dimensional Diffusive Growth 
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We consider a model ol aggregation, both diffusion-limited and ballistic, based on the Cayley tree. 
Growth is from the leaves of the tree towards the root, leading to non-trivial screening and branch 
competition effects. The model exhibits a phase transition between ballistic and diffusion-controlled 
growth, with non-trivial corrections to cluster size at the critical point. Even in the ballistic regime, 
cluster scaling is controlled by extremal statistics due to the branching structure of the Cayley tree; 
it is the extremal nature of the fluctuations that enables us to solve the model. 



I. INTRODUCTION 

Diffusion-limited aggregation (DLA) Q is a model of 
central importance in the field of fractal growth. For 
the theorist, it presents a model with very simple rules 
that gives rise to complex, branching patterns; for the 
experimentalist, it provides a model relevant for dielectric 
breakdown Q , viscous fingering || , electrodeposition || , 
and dendritic growth ||, among other physical processes. 

Recent progress in this field has been primarily re- 
stricted to d = 2 dimensions, due to the availability of 
powerful conformal techniques in the plane ||. How- 
ever, it remains important to elucidate the basic physics 
of DLA, branch competition and screening, for all di- 
mensions d. A phcnomenological model of branch com- 
petition was put forth by Halsey and coworkers 0. In 
this study, we present a high-dimensional model in which 
many of the branched growth ideas of Halsey may become 
exact. We find that the branched growth recursion equa- 
tions are controlled by extremal fluctuations. Finally, we 
comment on whether the extremal fluctuations are im- 
portant in DLA as well as in our model. 

For a fractal cluster of dimension D, the number of 
particles N(k) on a length scale k obeys 



N{k) oc k 1 



(1) 



It is known that for DLA, D > d — 1 g; it is believed 
that L> — > d — 1 as d — > oo. In order to study this high- 
dimensional limit, we choose to study diffusive and bal- 
listic growth on a Cayley tree. Unlike models of diffusive 
growth on Cayley trees in which the cluster grows from 
the root towards the leaves ||, our clusters grow from 
the leaves towards the root on a Cayley tree with k gen- 
erations of 2 k — 1 sites. The ballistic limit of this model 
has been considered previously There, the scaling 

behavior of N(k) was found only for a modification of 
the model, for which N(k) behaves differently. 

We observe non-trivial competition between the grow- 
ing branches, due to both geometrical and diffusive 
screening effects. We find that extremal fluctuations 
dominate the growth, leading to N(k) much less than the 
compact result of 2 k . N(k) still behaves exponentially in 
k, so that D is infinite. By examining the behavior of 



N(k) in different regimes, we find power law corrections 
to the exponential behavior, which may be interpreted as 
codimensions. 



II. THE MODEL 

Consider a Cayley tree with a constant branching fac- 
tor of 2 and a specified number of levels, k. We choose 
to initially occupy all the leaves of the tree, leaving all 
the other nodes empty. Therefore, a tree of k levels will 
have a total of 2 k — 1 nodes, of which 2 fe_1 — 1 are ini- 
tially empty. Random walkers will be released sequen- 
tially from the root of the tree to perform a random walk 
on the tree. When a random walker moves onto an oc- 
cupied site, it is added to the cluster at the site which it 
occupied previous to the last random walk step. 

Since each node has two daughters and one parent, 
there is an entropic bias to the random walk, causing it 
to drift towards the leaves. We offset this bias by defining 
a parameter 7 such that the walker moves towards a given 
daughter node with probability 1/(2 + 7), an d moves to- 
wards the parent node with probability 7/(2 + 7). At 
7 = 2, the walker is twice as likely to move to the parent 
as to any given daughter node, so the total probability 
of moving towards the leaves is equal to the probability 
of moving towards the root. At 7 — 0, the walker only 
moves towards the daughter nodes, so the growth is bal- 
listic, not diffusive. Growing branches can screen each 
other geometrically for all 7, but diffusive screening is 
present only for 7 > (see figure 1). 

Let us define the "absorbency" of a node. This is the 
probability that a random walker starting on that node 
will be absorbed into the cluster without first moving to 
that node's parent. The advantage of the absorbency is 
that it hides some of the non-local effects in the calcula- 
tion of random walk probabilities, and is useful for both 
analytical and numerical calculations. By definition the 
absorbency of an occupied node is equal to 1. 

The absorbency, a, of a node can be obtained solely 
from the absorbencies a, , a r of its left and right daughter 
nodes. We have 



7 + a, + a r 



(2) 
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For 7 < 2, we find that for a tree in the initial configura- 
tion, the absorbency of the root tends to a constant for 
large k. For 7 = 2 the initial absorbency of the root is 
equal to 1/fc. For 7 > 2, the absorbency of the root scales 
as (2/7) fe . This behavior of the absorbency suggests that 
7 = 2 is a critical point; wc confirm this suspicion below. 

Given the absorbencies, we can simulate the random 
walk as follows: start with the random walker at the 
root. Move the random walker towards the leaves at 
each step, until it hits an occupied node, with the relative 
probability of moving to the left or right daughters given 
by the relative absorbencies. 



III. BEHAVIOR AT 7 = 

We first consider the 7 = model, in which the random 
walker only moves towards the leaves, with equal proba- 
bility of moving left and right at each stage. Due to the 
branching nature of the Cayley tree, we can construct 
recursion relations to describe the growth. (These equa- 
tions will be very similar to those obtained phenomeno- 
logically in the branched growth model.) 

Define po(N,k) to be the probability that the root is 
unoccupied after N particles are added to a tree of k 
levels. Define pi(N,k) to be the probability that the 
root of the tree is occupied by the addition of the N-th 
particle. Define p2(N,k) to be the probability that the 
root was occupied before the iV-th particle was added 
(note that this means that some of the last particles must 
overflow from the tree) . Clearly, po (N, k) + p\ (N, k) + 
p 2 (N,k) = l. 

These probabilities obey simple recursion relations. 
Let N particles be added to a tree of j + 1 levels, and as- 
sume that the root of the tree is unoccupied at this point. 
All the particles lie below the root in the two subtrees of 
the tree, and with probability 2 _Ar (^,) we find that N' 
particles were added to the left subtree and N — N' were 
added to the right subtree. The root of a given subtree 
is either empty, or became occupied only as a result of 
the last particle to be added to that subtree. Therefore, 

N 



AT'=0 ^ ' 

(l-p 2 (N',j))(l-p 2 (N-N',j)). (3) 
Similarly, we can derive an equation for p\ . We find 

N 



Pi 



(N,j + 1)= 



N'=0 



2( Pl (N'-l,j))(l-p 2 (N-N',j)). 



(4) 



Combining the equations for p and p%, we obtain an 
equation for p 2 ■ 



Let us make a naive guess as to the solution of these 
equations. The guess, while incorrect, will be instruc- 
tive. The number of particles which go into each sub- 
tree will fluctuate. If we imagine that for N large these 
fluctuations become negligible, and further assume that 
Pi « po, we find that 



Po(N,j + l)=p (N/2,j) 2 



(5) 



The solution to this equation is po(N,j) — f(N/2 : >) such 
that f{x) = 1 for x < x c and f(x) = for x > x c . This 
means that the root of the tree is always empty if less 
than x c 2 k particles are added and after this the root is 
always occupied. 

Indeed, for j large, the noise in the recursion relation 
is exponentially small, so eq. (||) is not so different from 
eq. (|). However, if we perturb f(x) to f(x) = 1 — e(x) 
for x < x c , we find that e grows exponentially under eq. 
(||). From this we can guess that extremal statistics will 
become important. Since each node has an exponentially 
large number of nodes below it on the tree, even exponen- 
tially rare events may become important. In fact, using 
an argument based on extremal statistics, we now show 
that the average number of particles that must be added 
to a tree of k levels in order to occupy the root is given 
by N(k) « 2^2"^^. 

Consider a tree of k levels. Let k = I + m, and as- 
sume I » 1, m >> 1. We can imagine the tree as being 
made up of 2 l subtrees, each of to levels. Let us add 
(to — 1)2' particles to the tree. At this point, each of 
the subtrees has received roughly m — 1 particles (later, 
we will consider the effects of fluctuations in this number 
of particles). It is impossible for a subtree to complete 
its growth with only m — 2 particles added, since the 
subtree is to levels deep, with the leaves initially seeded. 
However, it is possible for one of these subtrees to com- 
plete its growth after only to — 1 particles are added. To 
compute this probability we note that the first particle 
added to the subtree can go anywhere; the second parti- 
cle must be added to the parent node of the first particle, 
which occurs with probability 2 - ( m ~ 2 ); the third parti- 
cle must go to the parent node of the second particle, 
with probability 2~( m ~ 3 ); and so on. This yields the re- 
sult that the filling probability, which is Pi(m — 1, to), is 

2— (m-l)(m-2)/2 ^ 2-»« 2 /2 

Thus, if 2'2~'" 2 / 2 is of order unity, it is likely that one 
of the subtrees will have been filled after only (to — 1)2' 
particles are added to the main tree. At this point, the 
next particle added from the root has probability 2~ l of 
attaching to the parent of the filled subtree, so that it 
will take approximately 2 l further particles added before 
the parent node of the subtree is occupied. At that point, 
it will take approximately 2 i_1 more particles before the 
grandparent of the subtree is occupied. Eventually, after 
2' + 2'~ 1 + 2'~ 2 + ... = 2 l+1 particles are added, we occupy 
the root of the whole tree. At this point, the total number 



2 



of particles added is (m — 1)2' + 2' +1 w m2 l . 

If we pick I and m correctly, we will thus take only m2 l 
particles to fill the tree. Since 2 l 2~ m I 2 must be of order 
unity, I w m 2 /2. For k >> 1, we find m w \/2ft and 
2 « fc — \/2fc. Once this number (m2 z ) particles has been 
added, we expect the tree to be filled so that 

N(k) < m2 l « 2 fc 2~ V2l v / 2£:. (6) 

Had we chosen 2, m such that 2 > m 2 /2 we would find 
many filled subtrees, and a weaker inequality than eq. 
. For I < m 2 /2, the argument does not work, because 

no subtree is filled. 

Now, suppose we add N(k) « 2 k 2~^ 2 ~ k \f2~k particles. 
If we pick 2 and m as above with m ~ v2fc, and if each 
of the subtrees receives of order N(k)/2 l « m parti- 
cles, then it is unlikely that any of the 2 l subtrees gets 
completed, hence the whole tree can not be completed. 
Thus if each subtree receives roughly N(k)/2 l particles, 
eq. (^) is an estimate and not just an inequality. 

Neglecting the effects of geometrical screening between 
subtrees, the number of particles arriving in a subtree fol- 
lows Poissonian statistics. Thus, while each subtree will 
receive roughly N(k)/2 l particles, there are rare events 
in which a subtree receives a large number of particles. 
Including these fluctuations, one can show that the prob- 
ability that a given subtree becomes occupied is 

oo n 

^ Pl (n,m)>e-"°, (7) 

n—l 

with no = N(k)/2 . The number of particles arriving in 
each subtree is still of order N(k)/2 l and the inclusion 
of these fluctuations does not change the estimate of the 
order of magnitude of N(k); the rare events considered 
in which a typical number of particles arrive in a sub- 
tree and form an atypical structure are much more likely 
than events in which an exponentially larger number of 
particles arrive and form a typical structure. 

To justify the neglect of geometrical screening, we note 
that for l,m, such that 2 < m 2 /2, we expect at most 
one subtree is occupied. Thus, the argument is self- 
consistent. 

Finally, let us express these estimates in terms of the 
recursion relations for Po,Pi- First, we have found that 
Po(N, m) = 1 for N < m — 1. However, 

p (m-l,m) = l-O(2- m2 / 2 ). (8) 

Next, we have used eq. @ to argue that po(2 l (m—l), m+ 
2) = l-0(2 z 2- m2 / 2 ) so that for 2'2~™ 2 / 2 w 1, the proba- 
bility po(2 (m — l),m+ 2) « 1, which indicates that the 
root is likely to be occupied. Using eq. (||) instead of eq. 
(||) corresponds to including the fluctuations expressed 
by eq. (@). 



IV. BEHAVIOR AT 7 > 

At 7 = 0, the only screening effect is the geometrical 
screening effect. We now repeat the argument based on 
extremal statistics for the case 7 > 0, when diffusive ef- 
fects begin to increase the competition between subtrees. 

For 7 > it is impossible to write closed recursion 
relations for the probabilities poj Pi ,P2- We can, however, 
write more general recursion relations, which include the 
absorbencies of the nodes. 

Still, for any 7, it is impossible to reach the root of 
an m level tree with m — 2 particles. The probability of 
reaching the root of an m level tree with m — 1 particles 
can be computed as in the 7 = case, now taking the 
absorbencies into account. For 7 < 2, the result is still 
given by eq. (|8j), while for 7 = 2, we find po((m— 1), m) = 
1 - 0(m! 2- m / 2 ), and for 7 > 2, P o((m - 1), m) = 1 - 
0(2~ m2/(21oS27 )). 

Using the same extremal argument as above, we find 
that, for 7 < 2, 

N(k) < 2 k 2- V2T V2k, (9) 

for 7 = 2, 

N(k) < 2 k 2-' /2 ~ k ~, (10) 

and for 7 > 2, 

N(k) < 2 fe 2-\ /2fel °S27 v / 2 fclog 2 7. (11) 

In deriving this, we assume that each of the 2 l subtrees 
continues to receive roughly N ' j2 l particles. However, 
now the number of particles arriving in a given subtree 
becomes non-Poissonian due to the effects of diffusive 
screening. Physically, this effect must reduce N(k). 

We still conjecture that eqs. (j^,[n)|,[ll]) are estimates, 
and not just inequalities. Suppose we add only N(k) << 
m2 l particles to the main tree, with m, I chosen to ob- 
tain the most stringent inequality. Let us show that this 
is not enough particles to grow to the root of the tree. We 
again do this by showing that it is justified to assume that 
each subtree receives roughly N(k)/2 l particles, so that 
none of the subtrees can be completely occupied. To jus- 
tify this assumption, we now need to show that diffusive, 
as well as geometrical, screening effects are negligible. 
When a subtree has received roughly N(k)/2 l << m par- 
ticles, the top 0(m) levels of each subtree are still likely 
to be completely empty. As a result, the absorbency 
of the subtree is is only slightly changed from that of a 
completely empty subtree and so the assumption of Pois- 
sonian statistics is reasonable. 

Interestingly, the 7 = 2 result for N(k) is reduced from 
the 7 < 2 result by an amount of order y/k, implying a 
codimension of 1/2 at the critical point. For 7 > 2, we 
find an infinite codimension, as the corrections to the 
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7 < 2 results are exponential, not power law. These 
results should be compared to the expected co-dimension 
of 1 for DLA in high dimensions. 



V. COMPARISON TO NUMERICS 

Using the absorbency technique, we have performed 
numerical calculations on the model for trees up to k = 24 
for 7 = 0, 1, 2, 4. In figure 2 we show log 2 (A r (fc)/2 fe ) aver- 
aged over 100 runs for each data point. As can be seen, 
N(k) << 2 k . Next, in figure 3, we show N(k) divided by 
the appropriate theoretical result for the scaling of N(k). 
For all curves, we see that the plot decreases initially, and 
then appears to stabilize at some given value, indicating 
that the theoretical result is in agreement with numerics. 
The theoretical calculation is only able to produce the 
asymptotic scaling of N(k). However, the prefactor for 
N(k) appears to be close to unity for 7 = 0. 



[5] J. Kertesz and T. Vicsek, J. Phys. A 19, L257 (1986). 

[6] M. B. Hastings and L. S. Levitov, Physica D 116, 244 
(1998); B. Davidovitch, H. G. E. Hentschel, Z. Olami, 
I Procaccia, L. M. Sander, and E. Somfai, Phys. Rev. E 
59, 1368 (1999); E. Somfai, L. M. Sander, and R. C. Ball, 
Phys. Rev. Lett. 83, 5523 (1999). 

[7] T. C. Halsey, Phys. Rev. Lett 72, 1228 (1994); T. C. 
Halsey and M. Leibig, Phys. Rev. A 46, 7793 (1992); T. 

C. Halsey, B. Duplantier, and K. Honda, Phys. Rev. Lett. 
78, 1719 (1997). 

[8] R. C. Ball and T. A. Witten, Phys. Rev. A 29, 2966 

(1984) . 

[9] J. Vannimenus, B. Nickel, and V. Hakim, Phys. Rev. B 
30, 391 (1984). 

[10] R. M. Bradley and P. N. Strenski, Phys. Rev. B 31, 4319 

(1985) . 

[11] C. Dasgupta and S. K. Ma, Phys. Rev. B 22, 1305 (1980); 

D. S. Fisher, Phys. Rev. Lett. 69, 534 (1992). 

[12] L. A. Turkevich and H. Scher, Phys. Rev. Lett. 55, 1026 
(1985). 



VI. CONCLUSION 

We have introduced a model that provides an inter- 
esting realization of some ideas of branched growth. We 
find in eq. (^) a simple recursion equation for the branch 
competition. However, the solution of this equation is 
unstable to noise. On the Cayley tree, the noise grows 
indefinitely leading to a situation controlled by extremal 
statistics, not unlike the "infinite disorder" critical points 
known in statistical mechanics JTT| . 

In finite dimensional DLA, we still expect that branch 
competition will "self-generate" noise, causing small per- 
turbations of a noiseless solution to grow, leading to a 
regime with non- vanishing fluctuations. Extremal statis- 
tics have been used to estimate the fractal dimension of 
DLA by Turkevich and Scher |]l2| . The relative impor- 
tance of large fluctuations as a function of dimensionality 
remains an unsolved problem. 
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FIG. 1. Screening for a tree with k — 4 after three growth 
steps. Filled circles represent occupied nodes. The empty 
circle represents a node that is geometrically screened with 
no growth possible on this site. For 7 > 0, there is a diffusive 
screening effect: the next random walker starting from the 
root is more likely to occupy the root as a result of a move 
onto the occupied left daughter node than it is to be added 
to the right half of the tree. 
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FIG. 2. log 2 (JV(fc)/2 fc )asafunctionofifc-lfor7 = 0,l,2,4 
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FIG. 3. N(k)/(2 k 2-' m V2k) as a 

function of fc - 1 fo r 7 = 0, 1; N{k)/{2 k 2- for 7 = 2; 
iV(fc)/(2' i: 2-\ /2fclo g27^2fclog 2 7) for 7 = 4. All curves asymp- 
tote to a constant. 
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